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We examine the dynamics of a skyrmion moving over a two-dimensional periodic substrate uti¬ 
lizing simulations of a particle-based skyrmion model. We specifically examine the role of the 
non-dissipative Magnus term on the driven motion and the resulting skyrmion velocity-force curves. 
In the overdamped limit, there is a depinning transition into a sliding state in which the skyrmion 
moves in the same direction as the external drive. When there is a finite Magnus component in 
the equation of motion, a skyrmion in the absence of a substrate moves at an angle with respect 
to the direction of the external driving force. When a periodic substrate is added, the direction 
of motion or Hall angle of the skyrmion is dependent on the amplitude of the external drive, only 
approaching the substrate-free limit for higher drives. Due to the underlying symmetry of the sub¬ 
strate the direction of skyrmion motion does not change continuously as a function of drive, but 
rather forms a series of discrete steps corresponding to integer or rational ratios of the velocity com¬ 
ponents perpendicular ((Vj_)) and parallel ((Vj|)) to the external drive direction: (Vj_)/(Vj|) = n/m, 
where n and m are integers. The skyrmion passes through a series of directional locking phases 
in which the motion is locked to certain symmetry directions of the substrate for fixed intervals of 
the drive amplitude. Within a given directionally locked phase, the Hall angle remains constant 
and the skyrmion moves in an orderly fashion through the sample. Signatures of the transitions 
into and out of these locked phases take the form of pronounced cusps in the skyrmion velocity 
versus force curves, as well as regions of negative differential mobility in which the net skyrmion 
velocity decreases with increasing external driving force. The number of steps in the transport 
curve increases when the relative strength of the Magnus term is increased. We also observe an 
overshoot phenomena in the directional locking, where the skyrmion motion can lock to a Hall angle 
greater than the clean limit value and then jump back to the lower value at higher drives. The 
skyrmion-substrate interactions can also produce a skyrmion acceleration effect in which, due to the 
non-dissipative dynamics, the skyrmion velocity exceeds the value expected to be produced by the 
external drive. We find that these effects are robust for different types of periodic substrates. Using 
a simple model for a skyrmion interacting with a single pinning site, we can capture the behavior 
of the change in the Hall angle with increasing external drive. When the skyrmion moves through 
the pinning site, its trajectory exhibits a side step phenomenon since the Magnus term induces a 
curvature in the skyrmion orbit. As the drive increases, this curvature is reduced and the side step 
effect is also reduced. Increasing the strength of the Magnus term reduces the range of impact 
parameters over which the skyrmion can be captured by a pinning site, which is one of the reasons 
that strong Magnus force effects reduce the pinning in skyrmion systems. 

PACS numbers: 75.47.Np,75.30.Kz,75.10.Hk,75.25.-j,75.75.-c 


I. INTRODUCTION 

Skyrmions were predicted to occur in certain magnetic 
systems! and were subsequently experimentally identified 
in the chiral magnet MnS^. Since this initial discovery 
there has been tremendous growth in the field as an in¬ 
creasing number of materials have been found that can 
support a skyrmion phased 0 . There are also numerous 
proposals on how to stabilize skyrmion states by utiliz¬ 
ing different materials properties or bilayers!!r— . Direct 
imaging of skyrmions with Lorentz microscopy! - —>!!! 
and other techniques ! 15 1 16 show that the skyrmions form 
a triangular lattice and have particle-like properties simi¬ 
lar to vortices in type-II superconductors!!. As an exter¬ 
nal magnetic field is increased, skyrmions emerge from a 
spiral state and their density initially increases and then 


decreases with field until the sample enters a uniform fer¬ 
romagnetic stateful. In bulk samples, skyrmions form 
three-dimensional (3D) line objects and occur in a limited 
range of fields and temperatures!!!, while for thin sam¬ 
ples the skyrmions exhibit two-dimensional (2D) proper¬ 
ties and are stable over a much larger range of fields and 
temperatures extending close to room temperature! - !. 
Skyrmions can be set into motion through the appli¬ 
cation of an external curren t 7 ' 10 : 18- !!, and it has been 
shown that there is a critical current above which the 
skyrmions depin into a sliding state!! - - 26 -. Skyrmion mo¬ 
tion can be directly observed with Lorentz microscopy!!! 
or deduced from changes in the transport properties, per¬ 
mitting the construction of effective skyrmion velocity 
versus applied force curves that show that the skyrmion 
velocity increases with increasing current 2 !. Other meth- 
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ods to move skyrmions include the use of temperature 
gradients^ - —, electric field s 32 ’ 33 , and coupling to a mag¬ 
netic tipi&. 

From an applications standpoint, skyrmions are at¬ 
tracting attention due to their potential use in racetrack 
memory devices where they would play a role similar to 
that of magnetic domain walls^i - —. Skyrmions have sev¬ 
eral advantages over domain walls due to their size and 
the fact that the current needed to depin a skyrmion 
can be orders of magnitude smaller than that needed 
to move domain wall s 20 ’ 34 . It has been experimentally 
demonstrated that individual skyrmions can be created 
or annihilated with a magnetic tip, indicating that it 
is feasible to read and write skyrmions^. Developing 
applications of skyrmions will require an understanding 
of how skyrmions interact with and move along tailored 
landscapes, so examining skyrmion dynamics on periodic 
substrates is an important step in this direction. 

Skyrmions have many similarities to vortices in type-II 
superconductors, such as the effective skyrmion-skyrmion 
interaction, which is repulsive and favors triangular or¬ 
dering, and the fact that both skyrmions and vortices 
can be driven by an external current. In the pres¬ 
ence of quenched disorder, skyrmions exhibit pinning- 
depinning transitions^, as observed in experiments^ 
and simulations^ - —, and these are similar to the 
depinning transitions observed for vortices in type-II 
superconductor s 17 ’ 37- 39 . There are important differences 
between the two systems, particularly the dominant role 
that non-dissipative effects can take in skyrmion motion. 
For superconducting vortices, non-dissipative effects such 
as a Magnus force are typically small, permitting the sys¬ 
tem to be effectively described as obeying overdamped 
dynamic s 17 ’ 38 . In a skyrmion system, by contrast, the 
Magnus term can strongly affect how the skyrmions in¬ 
teract with pinning sites and how they move under an 
external driv o 19 ’ 21 ’ 22 . Numerical simulations using con¬ 
tinuum and particle-based models for skyrmions inter¬ 
acting with pinning have shown that the Magnus term 
reduces the effective pinning in the system by creating a 
velocity component that is perpendicular to the force in¬ 
duced by a pinning site^^ -23 ’ 25 . As a result, skyrmions 
have a tendency to swing around the edge of a pinning 
site and escape, while for overdamped systems a particle 
moves toward the center of a pinning site and is much 
more likely to be pinned. 

Relatively little is known about how particles with a 
strong Magnus term move over a periodic substrate, and 
skyrmions are an ideal system to study such effects. In 
certain limits, a skyrmion can be effectively modeled as 
a point-like particle utilizing an equation of motion22*2iI 
derived from Thiele’s equation^. Comparison between 
continuum-based models and particle-based models of 
skyrmions moving in the presence of pinning have shown 
good agreemeni>22iS^. In this work we examine the dy¬ 
namics of a skyrmion moving over a square periodic sub¬ 
strate using a particle-based description given in Section 
II. We specifically examine the effect of changing the im¬ 


portance of the Magnus term relative to that of the damp¬ 
ing term. Despite the apparent simplicity of this system, 
we show in Section III that the Magnus term can induce 
a remarkably rich variety of dynamical behaviors that are 
absent in the overdamped limit. We find that the Hall 
angle for the skyrmion motion is dependent on the ex¬ 
ternal drive amplitude and approaches the substrate-free 
limit only at higher drives. Since the skyrmion is mov¬ 
ing over a periodic substrate, as the Hall angle changes 
the motion becomes locked to specific symmetry direc¬ 
tions of the substrate, producing a series of steps in the 
transport curves corresponding to integer and rational 
fractional ratios of the skyrmion velocity in the direc¬ 
tions parallel and perpendicular to the drive direction. 
At the transitions into the different directional locking 
phases, the skyrmion velocity exhibits a pronounced cusp 
accompanied by negative differential mobility in which 
the skyrmion velocity decreases with increasing external 
driving force. We map the extent of the locking phase as 
a function of the external drive and the ratio of the Mag¬ 
nus and dissipative terms, and find a rich structure of 
integer and fractional locking effects. We also describe a 
speedup effect for skyrmions interacting with a substrate 
where the skyrmion velocity can be higher than the ex¬ 
ternal drive. This effect is most prominent just above 
depinning and is caused by the Magnus term; it is absent 
in the overdamped case. In Section IV we consider a 
skyrmion scattering from a single pinning site for varied 
impact parameters to show how the Hall angle is reduced 
by the pinning due to a side-step phenomenon where the 
skyrmion trajectory is shifted by the Magnus term as the 
skyrmion moves through the pinning site. For increas¬ 
ing drive amplitude, the skyrmion trajectories become 
less curved and the size of the side step is reduced, so 
that the Hall angle approaches the clean limit for higher 
drives. Our results for the speedup effect and side steps 
are in agreement with recent theoretical and computa¬ 
tional studies by Miiller and Rosch, who considered a 
single skyrmion interacting with a defect site^S. In that 
work, the pinning potential is of a different form than the 
pinning sites we consider; however, the consistency of the 
two studies indicates that the Hall angle dependence on 
external drive and the speedup phenomenon are generic 
features of skyrmions interacting with pinning. 

There are other examples of particles moving over or¬ 
dered substrates, such as vortices in type-II superconduc¬ 
tors with periodic pinning arrays^ - — or colloids placed 
on optically created periodic substrate s 44 ’ 58 . In these sys¬ 
tems the dynamics is overdamped; however, there can be 
directional locking effects in which the particles preferen¬ 
tially move along symmetry directions of the underlying 
substrate as the direction of drive is rotated with respect 
to the substrate lattice^ - —. Such directional locking ef¬ 
fects can be exploited to perform particle separation in 
colloidal systems by setting up a system in which one 
particle species locks to the substrate while the second 
species does not, causing the two species to move at an 
angle with each othe r 47 ’ 48 . The direction of motion of 
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the locked particles undergoes a series of steps as a func¬ 
tion of the effective angle of drive with respect to the 
substrate. The steps are centered at integer and ratio¬ 
nal ratios of the angle of drive and form a devil’s stair¬ 
case structure ^ 48 ! 50 ! 51 ) 53 ! 54 . In the skyrmion system, we 
observe directional locking effects when the direction of 
external drive is fixed with respect to the substrate and 
the drive amplitude is varied. The directional locking and 
transitions into the locking states in the skyrmion system 
exhibit a number of features that have not been observed 
in overdamped systems, such as an overshoot effect in the 
locking direction and negative differential mobility at the 
locking transition. The steps in the transport response 
are also distinct from the Shapiro steps found in systems 
that can be effectively described as a particle moving 
over a periodic substrate under superimposed ac and dc 
drive s 55 ’ 56 . In the system we consider in this work, there 
is no imposed ac drive. 

II. SIMULATION AND SYSTEM 

We consider a skyrmion moving over a 2D square peri¬ 
odic substrate and utilize a particle-based description of 
the skyrmion from a recently developed equation of mo¬ 
tion for skyrmion s 22 ’ 23 . The dynamics of the skyrmion is 
obtained by integrating the equation of motion: 

r/R 

a d-^ = F i + «mZ X Vi + Ff . (1) 

Here = v, is the skyrmion velocity, and ad is a damp¬ 
ing term representing spin precession and dissipation of 
electrons localized in the skyrmion. 

The substrate force Ff is applied by placing N 
skyrmions in a square lattice with lattice constant a and 
fixing them in place. These skyrmions interact repul¬ 
sively with the mobile skyrmion. Such a substrate could 
be created experimentally by placing an array of mag¬ 
netic dots on the sample to create a background of ef¬ 
fectively immobile skyrmions. Alternatively, in recent 
continuum-based simulations it was shown that a de¬ 
fect site formed by removing a single spin creates a po¬ 
tential with a long range repulsion and a short range 
attraction^. Our system can then be viewed as con¬ 
taining a periodic array of such pinning sites in the 
limit where the moving skyrmion does not experience the 
shorter range attraction of the pins. We also consider a 
model in which the substrate potential is represented by 
a 2D sinusoidal array and find behavior very similar to 
that of a skyrmion moving through a pinned skyrmion 
lattice, indicating that our results are robust and cap¬ 
ture the general features of skyrmions moving over pe¬ 
riodic substrates. The force from the pinned skyrmions 
has the form Ff = where Rij is the 

distance between the driven skyrmion i and an immo¬ 
bile skyrmion j. Here K\ is the modified Bessel func¬ 
tion, which falls off exponentially for large Rij , and £ is a 
screening length which we take to be 1.0 in dimensionless 


units. The sample size is Lx L with L = 36, and the lat¬ 
tice constant of the substrate is a = 3.26. In the second 
system we model the substrate with the 2D periodic form 
F* = F p [cos 2 (7ra:/a)x + cos 2 (7r?//a)y], where a = 3.26 is 
the substrate lattice constant and F p is the amplitude of 
the substrate force. In both systems we employ periodic 
boundary conditions. 

The Magnus term F M = a m z x v,- produces a force 
that is perpendicular to the skyrmion velocity, where a m 
is the magnitude of the Magnus term. The driving force 
is Ff = Ff> d where d is the direction of the applied 
drive. Such a force could arise due to the application 
of an external current to the skyrmio n 22 ’ 23 . In most of 
this work, we take d = x. We measure the skyrmion 
velocity parallel, (Vj|), and perpendicular, (Vj_), to the 
drive. In the absence of a substrate or in the overdamped 
limit a m /ad = 0.0, the skyrmion moves in the direc¬ 
tion of the drive, while for a finite a m /ad the skyrmion 
moves at an angle 0 with respect to the drive, where 
0 = arctan((Vj_)/(Vji)) = arctan(a m /ad). Increasing 
otm/otd produces a larger angle 0 for the skyrmion mo¬ 
tion with respect to the external drive. To quantify the 
direction of motion we measure R = (Vj_)/(Vj|) so that 
the Hall angle 0 is given by arctan(I?). We increase the 
external drive in small increments of SFjj = 0.001 to 
0.005 and wait several thousand simulation time steps 
before measuring the average velocity in order to en¬ 
sure that the skyrmion velocity is in a steady state. We 
find that for smaller increments SFjj our results do not 
change. Throughout this work we impose the constraint 
a d + o%n = 1 in order to maintain a constant magnitude 
of the skyrmion velocity for varied ratios of a m /ad- 


III. VELOCITY-FORCE CURVES AND 
DIRECTIONAL LOCKING 

In Fig. |T] we plot (Vj|) and (Vj_) versus Fjj for a 
skyrmion moving over a periodic substrate in the over¬ 
damped case where ad = 1.0 and a m = 0.0. Figure 
[2] shows the system geometry, highlighting the motion 
of the mobile skyrmion through the background of po¬ 
tential maxima. In Fig. [T] there is a depinning tran¬ 
sition to a sliding state for Fd > 0.1 as indicated by 
(Vji) > 0. Here the skyrmion moves strictly in the di¬ 
rection of the applied drive so that (Vj_) = 0 for all Fjj 
and R = (Vj_)/(Vj|) = 0. The dashed line in Fig. |T] is 
the expected value of (Vj|) in the clean limit, and we find 
that as Fd increases, (Vj|) gradually approaches the clean 
limit value. Figure [2] shows the skyrmion trajectory for 
fixed Fjj = 0.5, where motion occurs in a straight line 
along the driving direction. 


A. Directional Locking 

In Fig. [3] we illustrate the transport behavior in a sam¬ 
ple with a m /ad = 0.45, a m = 0.41, and ad = 0.912085. 
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FIG. 1: (a) The velocity for a skyrmion under a dc driv¬ 

ing force Fd in the overdamped limit with oim/otd = 0 mov¬ 
ing over a periodic substrate. The drive is applied in the x- 
direction; the system geometry is illustrated in Fig. [2] (Vj|) is 
the velocity component parallel to the applied drive and (Vj_) 
is the velocity component perpendicular to the drive. Here 
there is a depinning transition into a sliding state where the 
skyrmion moves strictly in the direction of the applied drive. 
The Hall term R = (Vj_)/(Vj|) = 0.0 in this case. Dashed line: 
(Vj|) response in the absence of a substrate. 


y 


X 

FIG. 2: The geometry for the system in Fig.[T]with am/ad = 
0 where the skyrmion (large red dot) is driven in the x- 
direction under an applied drive Fd- The black line is the 
skyrmion trajectory and the smaller blue dots indicate the 
locations of the potential maxima in the periodic square sub¬ 
strate lattice. 
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FIG. 3: (a) (Vj_) and (Vj|) vs Fd for the same system as 
in Fig. [T] but with a m /ad = 0.45. (b) R vs Fd in the 

same sample. The dashed line indicates that in the clean 
limit R = —0.45. There are a series of dips in the velocity- 
force curves correlated with jumps in R, indicating that the 
skyrmion undergoes a series of transitions between locking to 
different symmetry directions of the substrate. The largest 
steps in R, corresponding to n/m ratios of 0/1, 1/4, 1/3, 2/5, 
3/7, and 1/2, are marked; there are also additional smaller 
higher order steps at other integer values of n and m. 




FIG. 4: The skyrmion trajectories and substrate maxima for 
the system in Fig. [3]with am/ad = 0.45. (a) |i?| = 0/1 state 
at Fd = 1.0. (b) \R\ = 1/4 state at Fd = 1.3. (c) |i?| = 2/5 
state at Fd = 2.2. (d) |i?| = 1/2 state at Fd = 2.75. 
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In the absence of a substrate the skyrmion would move 
at an angle of |0| = 24.227° with respect to the drive 
direction and would have R = (Vj_)/(Vj|) = —0.45. Fig- 
ureEIa) shows (Vj|) and (Vj_) vs Fd, and in Fig. [3^b) we 
plot R vs Ffj. Here there is a transition from a pinned to 
sliding state just above Fd = 0.1. The skyrmion moves 
strictly in the direction of the drive for 0.1 < Fd < 1.28, 
which corresponds to the R = 0 regime shown in Fig. EJb) 
for the same interval of Fd- Figure H|a) shows the 
skyrmion trajectory in the R = 0 region for Fd = 1.0 for 
the system in Fig. (Sja). The skyrmion no longer passes 
through the centers of the potential minima but its tra¬ 
jectory shifts closer to the bottom row of potential max¬ 
ima and develops an oscillation in the y direction that 
was absent in the a m /ad = 0 case shown in Fig. [2] As 
Fd is further increased the skyrmion shifts even closer to 
the row of potential maxima and its velocity parallel to 
the drive decreases, as indicated by the dip in (Vj|) near 
Fd = 1-28. As Fd increases further, both (Vj_) and (Vj|) 
increase, indicating that the skyrmion is now translating 
in the y direction as well as in the x direction. Additional 
dips in both (Vf|) and (Vj_) occur at higher Fd, with occa¬ 
sional regions containing multiple closely spaced smaller 
dips. Figure [3] indicates that these dips correlate with 
the jumps in R 1 and that between the dips, R remains 
constant. This means that the direction of motion or 
Hall angle of the skyrmion changes in a discrete fashion 
with increasing Fd ■ The steps in R appear at values that 
are rational ratios of (Vj_)/(Vj|) of the form n/m, where 
n and m are integers. In Fig. [3Jb) we highlight the steps 
at |i?| = 1/4, 1/3, 2/5, 3/7, and 1/2. There are also nu¬ 
merous additional steps in R for smaller intervals of Fd 
for higher order rational ratios of R = n/m. In general, 
the larger the values of n and in, the smaller the interval 
in Fd over which the step appears. 

During each step interval in R, the skyrmion follows 
an ordered periodic orbit, translating n substrate pla- 
quettes in the direction perpendicular to the drive for 
every m plaquettes it translates in the direction parallel 
to the drive. Figure IHa) shows the orbit of a skyrmion 
in the \R\ = 0/1 state, while Fig. Q|b) illustrates the 
orbit in the |i?| = 1/4 state at Fd = 1.3. Here the 
skyrmion moves periodically through the system at an 
angle 0 = arctan 1/4 = 14.036° with respect to the drive 
direction. In each period of the motion, the skyrmion 
translates to the right by four plaquettes in the drive 
or x direction and down by one plaquette in the per¬ 
pendicular or y direction, giving n/m = 1/4. At the 
transition out of the |i?| =1/4 state, the skyrmion slows 
down, producing a cusp in the velocity curves in Fig. Efb) 
and a dip in the net velocity near Fd = 1.327. A sim¬ 
ilar ordered orbit occurs in the |i?| = 1/3 state, where 
the skyrmion moves three plaquettes in the ^-direction 
and one plaquette in the y direction during every period. 
There is another dip in the net skyrmion velocity near 
Fd = 2.0 that occurs when the system transitions to 
the |i?| = 2/5 state. Figure Q|c) illustrates the orbit at 
this value of R for Fd = 2.2. Here the periodic orbit 



FIG. 5: A detailed view of the R. vs Fd curve for the same 
system in Fig. [3] at a m /otd = 0.45, where we examine the 
locking for drives up to Fd = 12.0. The system is locked to 
the \R\ = 1/2 value but then jumps to \R\ = 6/13 and 5/11 
before reaching |i?| = 4/9. As Fd is further increased, R 
gradually approaches the clean limit value, indicated by the 
dashed line. 


is more extended and the skyrmion moves two plaque¬ 
ttes in the negative y-direction for every five plaquettes 
in the ^-direction. Above the |i?| = 2/5 state, the sys¬ 
tem enters the |i?| = 3/7 state, followed by some higher 
order steps which occur over small intervals of Fd as 
shown in Fig. [5Ja) near Fd = 2.55. The system then 
reaches the |i?| = 1/2 state illustrated in Fig. |4|d) at 
Fd = 2.75, where the skyrmion moves two plaquettes in 
the ^-direction for every one plaquette in the y-direction. 
At the other higher-order locking steps, similar ordered 
orbits appear. 

The dips in (Vj|) and (Vj_) in Fig. [3] are associated 
with transitions in the skyrmion orbit from one direc¬ 
tionally locked state to another, with a corresponding 
change in the Hall angle. The symmetry of the square 
lattice determines the specific directions along which the 
skyrmion motion locks, so that for other geometries such 
as a triangular substrate, the locking directions will be 
different. At the transitions between locking steps, the 
net skyrmion velocity (V) = ((Vj|) 2 + (Vj,) 2 ) 1 / 2 decreases 
with increasing external drive Fd■ This phenomenon is 
known as negative differential mobility, and it has been 
observed in other systems where particles are driven over 
a periodic substrate, such as superconducting vortices 
moving over periodic pinning arrays where there are tran¬ 
sitions between different dynamical phases^. Negative 
differential mobility is also a common feature in semi¬ 
conductor devices^ and can be useful for creating logic 
devices. The ability to control differential mobility in 
skyrmion systems could open new approaches for appli¬ 
cations. 
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FIG. 6: (a,b) (Vj_) and (Vj|) vs Fd for (a) am/ad = 1-28 and 
(b) 1.91. (c,d) The corresponding R vs Fd curves. Dashed 
line in (d): the clean limit value for R. In (c) we highlight 
the \R.\ = 0/1 and |A| = 1/1 steps, and in (d) we highlight 
the \R\ = 1/1, 4/3, 3/2, 5/3, 7/4, 5/9 and 33/18 steps. 




FIG. 7: The skyrmion trajectories and substrate maxima, 
(a) The \R\ = 1/1 step for am/ad = 1-28 for the system in 
Fig- E|a,c). (b) The \R\ = 5/3 step for a m /ctd = 1.91 for the 
system in Fig. [6[b,d). (c) The |i?.| = 2/1 state for a m /ad = 
4.925 for the system in Fig. ifa,c). (d) The |i?| = 3/1 state 
B. Overshoot Effect at am/a d = 0.45 for arn / ad = 4.925 in Fig. HKc). 


At Fd = 3.0 in Fig. [3^b) , \R\ = 0.5, indicating the 
skyrmion is moving at a Hall angle of 0 = 25.565°, 
which is higher than the clean limit value of |i?| = 0.45 
or 0 = 24.23°. This phenomenon occurs when the 
clean value of 0 is oriented close to but slightly below a 
strong symmetry locking direction of the substrate. The 
skyrmion locks to the substrate and its Hall angle slightly 
exceeds that expected for the clean limit. As Fd is in¬ 
creased, the effectiveness of the substrate is reduced and 
the direction of motion of the skyrmion gradually ap¬ 
proaches the clean limit value, as shown in Fig. [5] where 
we plot a portion of the R vs Fd curve for the same sys¬ 
tem from Fig. EJb) but for drives up to a higher value 
of Fd = 12.0. The dashed line in Fig. [5] indicates the 
clean value limit of |i?| = 0.45. For Fd > 4.2, there is 
a jump from the \R\ = 0.5 = 1/2 state to a lower value 
of \R\ = 0.46143 = 6/13, followed by another jump to 
the \R\ = 0.4545 = 5/11 state. There is then a small 
region where the system locks to the \R\ = 0.45 = 9/20 
state before jumping to \R\ = 0.444 = 4/9. As Fd in¬ 
creases further, R gradually approaches the clean limit 
value, as indicated by the dashed line in Fig. [5jb). This 
shows that there can be an overshoot effect in the lock¬ 
ing behavior for a certain range of drives, giving rise to 
non-monotonic behavior of R as a function of Fd- We 
have observed similar overshoot effects for other values 
of / Q-d- 



FIG. 8: (a) (Vj|) and (V±) vs Fd at a m /aj = 4.925. (b) 

(Vj|) vs Fd for am/ad = 9.962. (c) R vs Fd for the sample in 
panel (a) with a m /ad = 4.925. The steps at \R\ = 2/1, 3/1, 
7/2, 4/1, 17/4, and 9/2 are highlighted, (d) R vs Fd for the 
sample in panel (b) with am/ad = 9.962, where the steps at 
\R\ = 5/1, 6/1, 7/1, 8/1, 17/2, and 9/1 are highlighted. 
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C. Higher Order Steps 


The specific features in the velocity force curves are 
strongly dependent on the value of a m /ad■ In Fig.|6|a,b) 
we show (Vj_) and (Vj|) for a m /ad = 1-28 and 1.91, re¬ 
spectively, while in Fig. [Gjc,d) we plot the correspond¬ 
ing R vs Fjj curves. For a m /ad = 1.28, the system is 
predominantly locked to the |i?| = 1/1 step as shown in 
Fig.[6jb). On this step, the skyrmion moves along the 45° 
direction as illustrated in Fig.[7fa) for F n = 1.0. Here the 
skyrmion follows a sinusoidal trajectory, moving by one 
plaquette in the ^-direction and one plaquette in the y- 
direction in a single period. At values of F D higher than 
those plotted in Fig. [S]A), further locking steps occur as 
|i?| approaches |i?| = a m /ad = 1.28. At a m /ad = 1.91, 
a larger number of steps occur, as indicated in Fig. Od) 
where we highlight the |i?| = 1/1, 4/3, 3/2, 5/3, 7/4, 
5/9, and 33/18 steps. For Fn values higher than those 
shown in Fig. [6jd) , additional steps appear as |i?| ap¬ 
proaches the clean limit value. Figure Ob) illustrates 
the skyrmion trajectory on the |i?| = 5/3 step for the 
system in Fig. Ob). The skyrmion moves five plaquettes 
in the direction perpendicular to the drive for every three 
plaquettes in the direction parallel to the drive during a 
single period. Fig. O a ) shows (Vj|) and (Vj_) versus Fn 
for a sample with a m /ad = 4.925, and the corresponding 
R vs Fn curve appears in Fig. O c )- I n this case, the 
clean limit Hall angle is large, and there is no longer a 
phase where the skyrmions move only in the direction of 
the drive, so the 0/1 step is lost. Instead, above depin- 
ning the motion jumps straight into the |i?| = 1.0 state. 
This is followed by steps at \R\ = 2/1, 3/1, 7/2, 4/1, 
17/4, and 9/2. There are also numerous smaller interme¬ 
diate steps corresponding to higher order fractions. In 
Fig.m c) we plot the skyrmion orbit on the |i?| = 2/1 
step for the system in Fig. HKc), while Fig. (TJd) shows 
the skyrmion orbit at |i?| = 3/1 for the same system. In 
general, as a m /ad increases, more steps become visible. 
In Fig.jHJb) we plot only (V||) versus Fn for a sample with 
ct m /&d = 9.962 to show more clearly the increased num¬ 
ber of features in the velocity-force curve. Figure [Hid) 
shows the corresponding R vs Fn curve, where we high¬ 
light the steps at |i?| = 5/1, 6/1, 7/1, 8/1, 17/2, and 
9/1. The largest observable integer step has a value of 
n that is the largest integer which is smaller or equal to 
the value of a m /ad- 


In Fig. [9])a) we plot the skyrmion trajectories at |i?| = 
5.0 for the system in Fig. (HJb,d) with a m /ad = 9.962, 
showing that the skyrmions circle around every fifth sub¬ 
strate maximum. At |i?| = 8.0, shown in Fig. EDb) , the 
trajectories are much straighter and the skyrmion moves 
at an angle that causes it to translate by eight plaquettes 
in the negative y direction for every one plaquette in the 
positive x-direction. 




FIG. 9: Skyrmion trajectories and substrate maxima for (a) 
the |A| = 5.0 state for the system in Fig. [Hlfod) at am/ad. = 
9.962, and (b) the |1?| = 8.0 state for the same system. 



FIG. 10: The regions of Fn vs am/ad where the pinned 
phase and the |i?| = 1/0, 1/1, 2/1, and 3/1 locking states 
occur. 


D. Arnol’d Tongues 

By conducting a series of simulations for varied a m /ad 
we can examine the evolution of the different locking 
phases. In Fig. [TO] we show the evolution of the pinned 
phase and the \R\ = 0/1,1/1, 2/1, and 3/1 locking phases 
for 0 < Fn < 3.0 and 0 < a m /ad < 5.0. The depinning 
threshold for the pinned phase remains roughly constant 
as a function of a m /ad■ For small a m /ad, the skyrmion 
motion is strictly in the x direction as indicated by the 
presence of the |i?| = 0/1 phase. As Fn increases, the 
range and the width of the lower order locking phases de¬ 
creases while new locking phases appear. The regions in 
which the locking occurs have the characteristic features 
of Arnol’d tongues, which occur in dynamical systems 
when there are two competing frequencie s 53 : 54 . In our 
system, the two frequencies are given by the inverse rate 
of translation along the x direction and the inverse rate of 
translation along the y direction when the system is on a 
locking step. These frequencies are quantized on the step 


































FIG. 11: (a) Fd vs am/ad for the a m /ad < 1.0 side of the 
\R\ = 1/1 locking step from Fig. 1101 (b) Fd vs a m /ad showing 
only the |i?| = 1/1 locking phase from Fig. [TO] The width of 
the step decreases for large Fd- 


due to the periodicity of the substrate; unlike the case of 
Shapiro steps, we do not apply any ac drive and the fre¬ 
quencies arise from the combination of dc motion and 
the substrate. If the substrate were rectangular instead 
of square, a different set of Arnol’d tongues would arise. 
At higher values of a m /ad , a larger number of locking 
steps appear in the velocity-force curves. For example, at 
otm/ctd = 2.5, as Fd increases the system passes through 
a small region of \R\ = 0/1 followed by the |i?| = 1/1 
locking phase and then by the |i?| = 2/1 phase. Another 
feature of the phase diagram is that at higher drives the 
width of each locking phase decreases. To illustrate this 
more clearly, we focus on the width of the \R\ = 1/1 step 
as a function of Fd, as highlighted in Fig. HHa) where 
we plot the location of the \R\ = 1/1 step over the nar¬ 
row range 0.8 < a m /ad < 0.96 out to Fd = 20, much 
higher than the maximum value of Fd shown in Fig. 1101 
The lower edge of the locking regime, which was dropping 
to lower values of a m /ad with increasing Fd in Fig. [10] 
bends back around for higher Fd as shown in Fig.QTj and 
shifts to higher a m /ad with increasing Fd- The resulting 
nose structure is the origin of the overshoot phenomenon, 
where over a certain range of a m /ad, the system locks 
to the |i?| = 1/1 step for lower Fd only to drop out of 
that step as Fd increases. Here the skyrmions can lock 
to the |i?| = 1/1 direction even though their resulting 
motion follows a higher angle 0 than would occur in a 
clean system. As Fd increases above the edge of the 
\R\ = 1/1 step, the skyrmion jumps to smaller values of 
\R\ in a series of steps. On each step the skyrmion locks 
to different symmetry directions which are closer to the 
clean system value of |i?| = a m /ad- The overall shape 
of the |i?| = 1/1 step is shown in Fig. [mb'), where we 
plot the range 0 < Fd < 25 and 0.8 < a m /ad < 1-32. 
There is a decrease in the extent of the locking phase at 
higher drives and at higher a m /ad values. The width of 
the |i?| =1/1 step gradually decreases and approaches a 
point centered at the a m /ad = 1.0 value for the highest 
drives. We find that the regions over which the other 
integer locking phases occur have similar shapes to that 
shown in Fig. ITTT bl. 



FIG. 12: (a) R vs Fd for the am/ad < 1.0 regime of the 
11?| = 1/1 locking phase from Fig. Illl bl. From top right 
to bottom right, am/ad = 0.8166, 0.8418, 0.8668, 0.98041, 
0.81475, and 0.9274. The |i?| = 1/1 step is marked, and the 
dashed line indicates the value of R for a m /ad = 0.8166 in the 
clean limit, (b) R vs Fd for the a m /ad > 1.0 regime of the 
|i?| = 1/1 locking phase from Fig. QjJ b). From top right to 
bottom right, a m /ad = 1.084, 1.1, 1.134, 1.17, 1.207, 1.246, 
and 1.2885. The \R\ = 1/1 step is marked and the dashed 
line indicates the value of R for a m/ad = 1.2885 in the clean 
limit. 


E. Overshoot Effect for Varied a m /ad 


A new set of locking steps arises for values of Fd 
greater than that where the |i?| = 1/1 locking occurs 
in Fig. HIT a). In Fig. IT2T al we plot R versus Fd for dif¬ 
ferent values of a m /ad on the a m /ad < 1.0 side of the 
\R\ = 1/1 locking phase shown in Fig. 11(b), where an 
overshoot effect occurs. The dashed line is the value of R 
at am/ad = 0.8166 in the clean limit; similar overshoots 
occur for the other values of a m /ad- The black curve 
shows that for a m /ad = 0.8166 in the presence of a sub¬ 
strate, the skyrmion locks to the |J?| = 1/1 direction for 
1.2 < Fd < 2.0 and then jumps to the \R\ = 4/5 = 0.8 
state. For Fd > 9.6, R jumps off of the \R\ = 0.8 step 
and, for higher values of Fd, further small jumps in R 
occur as R approaches the clean value limit. As a m /ad is 
increased toward a m /ad = 1.0, the width of the interval 
of Fd over which the system is locked to the \R\ = 1/1 
step increases since the clean limit Hall angle is closer to 
the \R\ = 1/1 angle of 45°. In each case, when the sys¬ 
tem jumps out of the \R\ = 1/1 step for increasing Fd, it 
can jump into a series of other locking steps as R grad¬ 
ually approaches the clean value. In Fig. [I2lbl we plot 
R versus Fd for the a m /ad > 1.0 side of the |i?| = 1/1 
locking phase from Fig.QTjb). The dashed line indicates 
the clean limit value of R for a m /ad = 1.2885. As a m /ad 
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FIG. 13: (a) The locking phases for Fd vs a m /ad highlight¬ 

ing the fractional locking steps with |J?| = 1/4, 1/3, 1/2, 2/3, 
3/4, 5/4, 4/3, 3/2, 5/3, and 7/4 steps along with the integer 
matching steps with \R\ = 0/1, 1/1, and 2/1. (b) Integer 
locking phases only for Fd vs a m /ad, for |1?.| =0/1 through 
|i?| = 12/1 from left to right. 


FIG. 14: (a) The net skyrmion velocity ( V} = ((Vj|) 2 + 

(Kl ) 2 ) 1 ' /2 vs Fd for am/ad = 9.962 (upper curve), 0.0 (lower 
curve), and the clean limit value (dashed line). Here there are 
regions of the am/ad = 9.962 curve where (V ) exceeds the 
clean limit value, indicating a speedup or acceleration effect, 
(b) AV = (V) - (H) c i ean for a m /a d = 9.962 (upper line), 
0.0 (lower line), and the clean limit value (dashed line). The 
speedup effect is indicated by regions in which AV > 0. 



approaches 1.0 from above, the extent of the region over 
which the system is locked to the |ii| = 1/1 state again 
increases. After jumping out of the |i?| = 1/1 state, the 
system jumps to larger values of \R\ in a series of smaller 
locking steps as it approaches the clean limit value of R. 
We observe similar changes in R for the other integer 
locking phases such as \R\ = 2/1 and |1?| = 3/1. 


F. Fractional Locking 

In Fig. 13(a) we plot the locking phases as a function 
of Fd and a m /ad over the range 0 < a m /ad < 2.5, high¬ 
lighting the fractional locking steps at the |1?| = 1/4, 
1/3, 1/2, 2/3, 3/4, 5/4, 4/3, 3/2, 5/3, and 7/4 states 
falling between the integer steps at \R\ = 0/1, 1/1, and 
2/1. Here the widths of the fractional locking states be¬ 
have similarly to those of the integer locking states. Fig¬ 
ure fl3T a) shows that for certain values of a m /ad, the 
fractional steps will be the dominant feature observed in 
transport. The higher order fractional steps, not shown 
in the figure, exhibit similar features. The overshoot ef¬ 
fect described for the \R\ = 1/1 locking step in Fig. [TT] 
occurs for all of the integer and fractional locking steps as 
well, with the overall width of each of the steps decreasing 
with increasing Fd- In Fig. 1131 b) we plot only the integer 
locking regions from |1?| = 0/1 to |1?| = 12/1 as a func¬ 
tion of Fd and a m /ad for the range 0 < a m /ad < 16.5. 
Between each of the integer steps, there is a series of 
fractional steps (not shown) similar to that illustrated in 
Fig-EHa). The fractional steps have the form N + n/m, 


where N are the integer steps. 

G. Speedup Effects 

The Magnus term can produce an acceleration or 
speedup effect of the skyrmion, in which the speed of 
the skyrmion is higher in the presence of a substrate 
than it would be in the absence of a substrate or in 
the overdamped limit. In Fig. [Uf a) we plot the net 
skyrmion velocity (F) = ((Vj|) 2 + (Vj_) 2 ) 1 / 2 versus Fd for 
dm/cud = 9.962, the overdamped case of a m /ad = 0, and 
the clean limit. In the overdamped case, ( V) is always 
smaller than the clean limit value, indicating that the 
substrate does not accelerate the skyrmion. In contrast, 
for a m /ad = 9.962 there are clear regions where ( V ) is 
higher than the clean value. This effect is most promi¬ 
nent for values of Fd just above the depinning thresh¬ 
old. To show the speedup more clearly, in Fig. iMT h) 
we plot AV = (V) — (F) c iean for a m /ad = 9.962 and 
0.0. For the overdamped case, AV < 0.0 over the en¬ 
tire range of Fd- The lowest values of AV occur just 
at depinning, and then AV gradually approaches zero as 
Fd increases. The behavior of the overdamped system 
is similar to the velocity-force curves observed in other 
overdamped systems such as colloids driven over periodic 
substrate array s 44 ' 58 . For a m /a.d = 9.962, AV shows a 
series of peaks which are correlated with transitions be¬ 
tween the different locking regimes. For 0.1 < Fd < 0.44, 
AH > 0.0 due to the speedup effect, where there is an en- 
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FIG. 15: (V) vs a m /ad for fixed Fd = 0.2. The dashed 

line indicates the clean limit value of (V) = 0.2. The speedup 
effect occurs when ( V) > 0.2 and increases in magnitude with 
increasing ctm/ctd- 


hancement in the net velocity of up to twice the velocity 
in the clean limit. As Fjj is further increased the mag¬ 
nitude of the speedup effect decreases and there are sev¬ 
eral intervals of Fq where the skyrmion is moving signifi¬ 
cantly slower than it would in the clean limit. For further 
increases in Fjj, AV approaches zero. The speedup effect 
is generated by the non-dissipative terms in the equation 
of motion which cause the skyrmions to be accelerated 
through certain portions of the substrate potential. 

In Fig. [15] we plot (V) versus a m /ad for a fixed Fd = 
0.2, where the dashed line indicates the clean limit value 
of (V) = 0.2. Values of (V) > 0.2 indicate a speedup 
effect. For a m /ad = 0, in the overdamped limit, (V) < 
0.2. We find a series of oscillations in ( V) corresponding 
to the different locking phases through which the system 
passes as a function of a m /a.d■ Here, the magnitude of 
the speedup effect increases on average with increasing 

A speedup effect for a driven skyrmion interacting with 
a defect has also been observed in simulations by Muller 
and Rosch, who find that for some cases the defect causes 
a net increase in the skyrmion velocity^. They also 
find that the magnitude of the speedup effect decreases 
when an externally imposed skyrmion drift velocity is in¬ 
creased. This is similar to what we observe for the peri¬ 
odic substrate case where for increasing external driving 
force the speedup effect is reduced as shown in Fig- fTH bb 
In Refi2& the speedup was up to an order of magnitude 
higher than the velocity in the clean limit, while we ob¬ 
serve a velocity enhancement of only up to a factor of 2. 
This is because in our system there is a minimum critical 
force required to depin the skyrmion, whereas in Refj^ 
there was no lower bound on the imposed drift velocity. 
Since the speedup is increased for lower drives, Muller 
and Rosch could access lower drives and obtain larger 
velocity enhancements from the speedup effect. 



FIG. 16: A gray scale image of a portion of the 2D analytic 
substrate. The external drive is applied along the ^-direction. 



FIG. 17: The transport for a skyrmion moving over the 2D 
sinusoidal substrate illustrated in Fig. [16] for F v = 1.5 and 
Om/ad = 4.925. (a) The net skyrmion velocity ( V) vs Fd 
shows a series of cusps. Dashed line: the clean value limit, 
(b) (Vx) and (Vj|) vs Fd for the same system, (c) R vs Fd for 
the same system with the |A| = 1/1, 2/1, 3/1, and 4/1 steps 
highlighted. There are also numerous smaller scale fractional 
locking steps. 


H. Two Dimensional Analytic Substrate 

In order to understand how general our results are for 
different details of the periodic substrate, we next con¬ 
sider a 2D analytic substrate with the same lattice con¬ 
stant a used in obtaining the previous results. The force 
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from the substrate is given by F s = F p [cos 2 (nx / a)x + 
cos 2 (iry /a)y], where F p is the maximum force exerted by 
the substrate and a = 3.26 is the lattice constant as in 
the system shown in Fig. [2] Figure [lB] shows a gray scale 
of a portion of the substrate, where the potential max¬ 
ima are highlighted. In Fig. m a) we plot (V) versus Fd 
for a skyrmion moving over a substrate with F p = 1.5 
for a m /ad = 4.925. The dashed line indicates the clean 
limit value of (V). Figure fTTl bl shows (Vj|) and (Vj_) vs 
Fjj for the same system, while in Fig. HUT) we plot the 
corresponding R vs Fp curve with a dashed line indicat¬ 
ing the value of I? in a clean system and with the locking 
steps at \R\ = 1/1, 2/1, 3/1, and 4/1 highlighted. Al¬ 
though there are some differences in the details, we find 
the same general transport features for the analytic sub¬ 
strate as we observed for the pinned skyrmion substrate 
at OL m lad = 4.925, including the locking of the skyrmion 
motion to different symmetry directions and the steps 
in R at integer and fractional ratios of (Vj_)/(Vj|). Fig¬ 
ure fTTT al shows that there are a similar series of dips 
in ( V) corresponding to transitions between the different 
locking phases. We find that with an analytic potential, 
a somewhat larger number of small looking steps can be 
resolved, as shown for 1.25 < Fp < 1.75 in Fig. flTT ab 
and the steps in R in Fig. E(c) for this region have a 
devil’s staircase structure. Figure [Tri al also shows that 
(V) exhibits regions where it is higher than the clean limit 
(dashed line), indicating that the same type of speedup 
effect occurs for the analytic potential. We have also ex¬ 
amined the locking effects on the analytic potential for 
other values of a m /ad and F p . and find that all the fea¬ 
tures highlighted in Fig. [17] are robust. This indicates 
that the directional locking effect is a generic feature of 
skyrmions moving over periodic substrates. 


IV. SCATTERING OFF A SINGLE PINNING 
SITE 

In order to better understand the dependence of the 
Hall angle on the external drive, we consider the case of 
a skyrmion scattering from a single pinning site. We 
drive the skyrmion toward the pinning site for varied 
impact parameters b, and measure the resulting shift in 
the skyrmion position perpendicular to the driving force 
for the outgoing state. We model the pinning site as a 
parabolic trap with radius R p and a maximum force of 
F p . The skyrmion is driven from a point outside the trap 
towards the trap with an external drive Fp applied in 
the positive y-direction. In Fig. IlSf al we show the over¬ 
damped case with a m /ad = 0.0, F p = 0.1, R p = 0.35, 
and Fp = 0.05, We highlight the trajectory of a skyrmion 
with an impact parameter of b = 0.0. For this value of 
Fp. all skyrmion trajectories that contact the pinning 
site form straight lines directed toward the equilibrium 
position of the pinned state along the x = 0.0 line. The 
pinned equilibrium position is shifted in the positive y- 
direction from the center of the pinning site due to the 



FIG. 18: Skyrmion trajectories (lines) for a single skyrmion 
at different impact parameters b interacting with a single pin¬ 
ning site of radius R p = 0.35 and maximum pinning force 
of F p for an external drive Fd applied in the positive y- 
direction. (a) am/ad = 0.0, F p = 0.1, and Fd = 0.05. The 
skyrmions that interact with the pinning site are captured, 
(b) am/ctd — 0, F p = 0.1, and Fd = 0.12. The skyrmions can 
escape from the pinning site, (c) am/ad = 10.0, F p = 0.1, and 
Fd = 0.05. The dashed line is a trajectory that a skyrmion 
would follow in the absence of the pinning site. For a cer¬ 
tain range of impact parameters, the skyrmion is captured. 
For other impact parameters, skyrmions that escape from the 
pinning site have their trajectories shifted as highlighted by 
the thick line, (d) am/ad = 10.0, F p = 0.1, and Fd = 0 .12. 
Here the shift in the trajectories of the skyrmions that pass 
through the pinning site is reduced and there is less curvature 
in the trajectories. 

applied external force. Figure fTST b') shows the same sys¬ 
tem with higher Fp = 0.12. The trajectories are shifted 
toward the center of the pinning site, but since Fp > F p 
the skyrmions can escape from the pinning site. The tra¬ 
jectory for an impact parameter 5 = 0 comes in and out of 
the pinning site at x = 0.0, so there is no shift, while the 
shift of the trajectories for non-zero impact parameters 
are symmetric across x = 0. 

Fig. 18(c) shows the trajectories for the same system 
at am/ad = 10.0 for Fp = 0.05, where the external drive 
is still applied in the positive y-direction. In this case, 
the trajectories that do not intersect with the pinning 
site move at an angle of | arctan(a m /ad)| = 84.3° with 
respect to the y-axis. The dashed line represents the 
trajectory a skyrmion would follow in the absence of the 
pinning site. When the skyrmion moves through the pin¬ 
ning site, the trajectory is no longer straight as in the 
overdamped case but is now strongly curved due to the 
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Magnus force. The skyrmions that encounter the pinning 
site on the upper left side become trapped by entering 
a spiral trajectory that carries them to the equilibrium 
pinned position in the upper portion of the pinning site. 
Skyrmions that enter the pinning site on the lower left 
side can escape from the pin; however, the position of the 
outgoing skyrmions are shifted in the positive y direction 
as indicated by the highlighted trajectory in Fig. 0c). 
This shift is similar to the side jump effect that occurs for 
electron scattering^, where the interaction with pinning 
or disorder shifts the outgoing trajectory of the particle 
relative to its incoming trajectory. In the case of a pe¬ 
riodic substrate, a skyrmion would be repeatedly shifted 
as it moves through the system, so that for a m /ad = 10 
the average direction of skyrmion motion follows a Hall 
angle that is less than the clean value limit of 84.3°. Fig¬ 
ure ante) also shows that when the skyrmion is captured 
it becomes pinned along the x = 0.0 line at a position 
above the center of the pin. The shift of the pinned equi¬ 
librium position to this location is caused by the bias 
from the external drive, and the location of the equi¬ 
librium position is not changed by the inclusion of a fi¬ 
nite Magnus term. Figure [TBI dl shows the trajectory for 
oim/ctd = 10 with a higher drive of Fp = 0.12. Here, the 
size of the side jump is reduced, so that the Hall angle 
that would be observed for motion through a periodic 
substrate is closer to the pin-free value of 84.3°. This in¬ 
dicates that the Hall angle increases with increasing driv¬ 
ing force, as also observed in the simulations with a peri¬ 
odic substrate. The shift is reduced at the higher drives 
in Fig. I18f di. since the skyrmion is moving more rapidly 
through the pinning site and the trajectory spends less 
time being bent. 

To quantify the the dependence of the magnitude Sr of 
the shift or side jump of the trajectory as a function of 
the external drive, in Fig. [TO] we plot Sr for a skyrmion 
approaching the pinning site with 6 = 0. Figure I19l al 
shows the shift for F p = 0.1, R p = 0.35, and a m /ad = 
10. The shift is highest at low Fp, starting near Sr = 
0.3 in Fig. HHa) and gradually approaching zero as Fp 
increases. In Fig.[l9|b) we plot Sr versus Fp for a sample 
with F p = 0.1, R p = 0.35, and a m /ad = 1-0. Here the 
initial value of Sr is smaller due to the smaller value of 
otm/otd ■ From an initial value of Sr = 0.15, the magnitude 
of the shift decreases with increasing drive. For a m /ad = 
0.0 (not shown), Sr = 0 for all values of Fp. If we apply 
a fit to the decrease of the shift as a function of drive, 
we find Sr oc F^ v , with v = 1.44 — 2.0 depending on 
the choice of the low drive cutoff. In the work of Muller 
and Rosch for a skyrmion scattering off a single defect^, 
they found through both simulations and perturbation 
that the shift as a function of drive goes as Sr oc Fd- 
Even though the pinning potential we use has a different 
form than that used in Ref^, our results are consistent 
with the pinning site inducing a trajectory shift. 

In order to understand the shift of the skyrmion and 
when the skyrmion will be captured by a pinning site, in 
Fig.[2D]we plot Sr as a function of the impact parameter 6. 



FIG. 19: Sr, the magnitude of the trajectory shift, vs Fd 
for a skyrmion entering a single pinning site with an impact 
factor of zero, (a) F p = 0.1, R p = 0.35, and am/ad = 10.0. Sr 
approaches zero at higher drives, (b) The same for am/ad = 
1 . 0 . 
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FIG. 20: The trajectory shift Sr vs impact parameter b for 
a skyrmion moving through a pinning site with F p = 0.1 and 
R p = 0.35. (a) am/ad = 0 for Fd = 0.12, 0.16, and 0.20, 
from upper left to lower left. Sr is symmetric across 6 = 0. 
(b) am/ad = 1.0 for F D = 0.085, 0.09, 0.10, 0.11, 0.12, 0.14, 
0.16, 0.18, and 0.20, from upper left to lower left. Here the 
shift is asymmetric across 6 = 0. For values of 6 where there 
are no points, the skyrmion is captured by the pinning site, 
(b) am/ad = 10 for Fd = 0.03, 0.05, 0.08, 0.12, and 0.20, 
from upper left to lower left. The shifts become larger with 
increasing a m/ad- 
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Figure EOT a) shows the overdamped case a m /ad = 0.0 for 
F p = 0.1 and r p = 0.35 at F D = 0.12, 0.16, and 0.2. At 
b = 0, Sr = 0 and there is no shift in the trajectory, while 
shifts for positive and negative values of b are symmetric 
across b = 0 so that the integrated shift over all impact 
parameters is zero. For Fn < F p = 0.1 the skyrmion 
is always captured by the pin. Figure I20f bl shows Sr 
versus b for a sample with a m /ad = 1.0 at drives ranging 
from Fn = 0.085 to Fn = 0.20. Here, Sr decreases with 
increasing Fn- and the shifts are asymmetric for positive 
and negative values of b, so that the integrated shift over 
all impact parameters is positive. At Fn = 0.085 and 
0.09, values of b at which there are no points on the curve 
indicate that the skyrmion was captured by the pinning 
site. For higher values of Fn , |<5r| is larger for b < 0 than 
for b > 0. The results in Fig.jSOJb) show that increasing 
otm/oib reduces the range of impact parameters where an 
incoming skyrmion is captured by the pin, and that for 
some ranges of b a skyrmion can escape the pin even when 
Fn < F p , in contrast to the overdamped case where the 
skyrmion is always trapped whenever Fn < F p . 

Figure EUT c'l shows Sr versus b for a m /ad = 10.0 over 
the range Fn = 0.03 to Fn = 0.2. Here, the ability of the 
pinning site to capture the skyrmion is further reduced, 
and the skyrmion does not become pinned for more than 
half the impact parameters at Fn = 0.03 even though 
this drive is substantially smaller than the maximum pin¬ 
ning force F p = 0.1. At Fn = 0.05 a skyrmion can only 
become captured if b > 0.25. The overall shifts are posi¬ 
tive over a much wider range of Fn, and the magnitude 
of Sr decreases with increasing Fn- This result shows 
that inclusion of a Magnus term in the dynamics reduces 
the ability of pinning sites to capture particles since the 
Magnus term induces a side jump or shift that permits 
the particles to escape from the pin. 


V. CONCLUSION 

We have numerically examined a skyrmion under a dc 
drive moving over a two-dimensional square pinning sub¬ 
strate for varied ratios of the Magnus force term a m to 
the damping ad- In the overdamped limit a m /ad = 0, 
there is a single depinning transition into a sliding state 
where the skyrmion moves in the direction of the applied 
drive. For a finite a m /ad we find that the skyrmion 
direction of motion or Hall angle depends on the magni¬ 
tude of the external drive and gradually approaches the 
substrate-free limit at high drives. Due to the symme¬ 
try of the underlying substrate, the Hall angle does not 
change continuously but passes through a series of steps 
as the skyrmion motion becomes locked to certain sym¬ 
metry directions of the substrate. These steps occur at 
integer and rational fractional ratios n/m of the perpen¬ 
dicular to parallel velocity components of the skyrmion 
motion with respect to the direction of the applied drive, 
where n and m are the number of plaquettes the skyrmion 
moves in the perpendicular and parallel directions, re¬ 


spectively. We find that the Hall angle generally increases 
with increasing external drive, but that there can be an 
overshoot effect in which the Hall angle is larger than ex¬ 
pected for the substrate-free or clean limit. In this case, 
as the drive increases the Hall angle jumps back to a lower 
value closer to the clean limit value. At the transitions 
between different directional looking steps, the skyrmion 
velocity shows a striking series of cusps or dips where 
the skyrmion slows down for increasing Fn, producing 
a negative differential mobility at the transitions. As 
o-m/otd increases, the number of transitions between dif¬ 
ferent looking steps increases, as we map out in a series of 
phase diagrams. The directional locking effects exhibited 
by the skyrmions are very distinct from the directional 
locking effects previously observed for overdamped parti¬ 
cles such as vortices and colloids interacting with periodic 
substrates. In the overdamped system, the angle of the 
external drive must be changed with respect to the sym¬ 
metry direction of the substrate in order to induce lock¬ 
ing steps, whereas for the skyrmion system the driving 
direction remains fixed; only the magnitude of the driving 
force is changed. Additionally, the overdamped systems 
exhibit neither the negative mobility phenomenon at the 
transitions between steps nor the overshoot effect. 

We find that the skyrmion motion can exhibit a 
speedup or acceleration effect where the interaction of the 
skyrmion with a pinning site can accelerate the skyrmion 
such that the skyrmion velocity is higher than the value 
that would be induced by the external drive alone. This 
effect is generally enhanced at the lower drives and sup¬ 
pressed near the directional locking transitions. It con¬ 
trasts strongly with the behavior of overdamped systems, 
where interactions with pinning or a substrate always de¬ 
crease the velocity of an overdamped particle. We find 
that all of these features are robust for different forms of 
the periodic substrate. 

To better understand how the pinning can induce a 
Hall angle dependence on the magnitude of the drive, we 
consider a skyrmion scattering from a single pinning site. 
In the overdamped limit with a m /ad = 0, the skyrmion 
becomes trapped by the pin whenever the external drive 
is less than the magnitude of the maximum pinning force 
that can be exerted by the pinning site. For external 
drives larger than this value, the skyrmion escapes from 
the pin and its trajectory is shifted perpendicular to the 
direction of skyrmion motion. The shift is symmetric for 
positive and negative impact parameters, so there is no 
net shift of the trajectory when all impact parameters 
are integrated. In contrast, at a finite a m /ad , the mag¬ 
nitude of the trajectory shift depends asymmetrically on 
the impact parameter, producing a nonzero net shift of 
the skyrmion as it moves through the pinning site when 
all impact parameters are integrated. This shift is similar 
to the side jump phenomenon found in electron scatter¬ 
ing. As the external drive is increased, the net shift de¬ 
creases and the skyrmion motion approaches the pin-free 
trajectory track. The trajectory shift is responsible for 
the modified Hall angle experienced by the skyrmion at 
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low drives, and the gradual decrease in shift with increas¬ 
ing drive causes the Hall angle to gradually approach the 
clean limit value. For the case of a periodic substrate, the 
additional symmetry of the substrate prevents the angle 
of the skyrmion motion from changing continuously; in¬ 
stead, the motion changes in a series of steps with in¬ 
creasing drive as it approaches the clean limit value. We 
also find that for finite a m /ad , depending on the im¬ 
pact parameter the skyrmion may or may not become 
pinned regardless of whether the maximum pinning force 
is smaller than the driving force, and that for increas¬ 
ing a m /ad , the range of impact parameters over which 
the particle is not pinned even when the driving force is 


smaller than the maximum pinning force increases. This 
is one of the reasons that pinning of skyrmions is much 
weaker than pinning in overdamped systems such as su¬ 
perconducting vortices. 


Acknowledgments 

This work was carried out under the auspices of the 
NNSA of the U.S. DoE at LANL under Contract No. 
DE-AC52-06NA25396. 


1 A. Bogdanov and A. Hubert, J. Magn. Magn. Mater. 138, 
255 (1994); U. K. Rofiler, A. N. Bogdanov, and C. Pflei- 
derer, Nature (London) 442, 797 (2006). 

2 S. Miihlbauer, B. Binz, F. Jonietz, C. Pfleiderer, A. Rosch, 
A. Neubauer, R. Georgii, and P. Boni, Science 323, 915 
(2009). 

3 X.Z. Yu, Y. Onose, N. Kanazawa, J.H. Park, J.H. Han, 
Y. Matsui, M. Nagaosa, and Y. Tokura, Nature 465, 901 
( 2010 ). 

4 X.Z. Yu, N. Kanazawa, Y. Onose, K. Kimoto, W.Z. Zhang, 
S. Ishiwata, Y. Matsui, and Y. Tokura, Nature Mater. 10, 
106 (2011). 

5 S. Seki, X.Z. Yu, S. Ishiwata, and Y. Tokura, Science 336, 
198 (2012). 

6 S.X. Huang and C.L. Chien, Phys. Rev. Lett. 108, 267201 

( 2012 ). 

7 X.Z. Yu, N. Kanazawa, W.Z. Zhang, T. Nagai, T. Hara, 
K. Kimoto, Y. Matsui, Y. Onose, and Y. Tokura, Nature 
Commun. 3, 988 (2012). 

8 S. Heinze, K. von Bergmann, M. Menzel, J. Brede, A. Ku- 
betzka, R. Wiesendanger, G. Bihlmayer, and S. Bltigel, 
Nature Phys. 7, 713 (2011). 

9 N. Nagaosa and Y. Tokura, Nature Nanotech. 8, 899 
(2013). 

10 X.Z. Yu, Y. Tokunaga, Y. Kaneko, W.Z. Zhang, K. Ki¬ 
moto, Y. Matsui, Y. Taguchi, and Y. Tokura, Nature Com¬ 
mun. 5, 3198 (2014). 

11 A.B. Butenko, A.A. Leonov, U.K. Rofiler, and A.N. Bog¬ 
danov, Phys. Rev. B 82, 052403 (2010). 

12 R.L. Silva, L.D. Secchin, W.A. Moura-Melo, A.R. Pereira, 
and R.L. Stamps, Phys. Rev. B 89, 054434 (2014). 

13 S. Polesya, S. Mankovsky, S. Bornemann, D. Kodderitzsch, 
J. Minar, and H. Ebert, Phys. Rev. B 89, 184414 (2014). 

14 S. Banerjee, J. Rowland, O. Erten, and M. Randeria, Phys. 
Rev. X 4, 031045 (2014). 

15 N. Ronning, C. Hanneken, M. Menzel, J.E. Bickel, B. 
Wolter, K. von Bergmann, A. Kubetzka, and R. Wiesen¬ 
danger, Science 341, 636 (2013). 

16 P. Milde, D. Kohler, J. Seidel, L.M. Eng, A. Bauer, A. 
Chacon, J. Kindervater, S. Miihlbauer, C. Pfleiderer, S. 
Buhrandt, C. Schiitte, and A. Rosch, Science 340, 1076 
(2013). 

17 G. Blatter, M.V. Feigelman, V.B. Geshkenbein, A.I. 
Larkin, and V.M. Vinokur, Rev. Mod. Phys. 66, 1125 
(1994). 


18 F. Jonietz, S. Miihlbauer, C. Pfleiderer, A. Neubauer, W. 
Miinzer, A. Bauer, T. Adams, R. Georgii, P. Boni, R.A. 
Duine, K. Everschor, M. Garst, and A. Rosch, Science 330, 
1648 (2010). 

19 J. Zang, M. Mostovoy, J.H. Han, and N. Nagaosa, Phys. 
Rev. Lett. 107, 136804 (2011). 

20 T. Schulz, R. Ritz, A. Bauer, M. Haider, M. Wagner, 
C. Franz, C. Pfleiderer, K. Everschor, M. Garst, and A. 
Rosch, Nature Phys. 8, 301 (2012). 

21 J. Iwasaki, M. Mochizuki, and N. Nagaosa, Nature Com¬ 
mun. 4, 1463 (2013). 

22 S.-Z. Lin, C. Reichhardt, C.D. Batista, and A. Saxena, 
Phys. Rev. B 87, 214419 (2013). 

23 S.-Z. Lin, C. Reichhardt, C.D. Batista, and A. Saxena, J. 
Appl. Phys. 115, 17D109 (2014). 

24 R.E. Troncoso and A.S. Nunez, Phys. Rev. B 89, 224403 
(2014). 

25 C. Reichhardt, D. Ray, and C.J. Olson Reichhardt, 
arXiv: 1409.5457 

26 J. Miiller and A. Rosch, arXiv:1411.2857 

27 K. Everschor, M. Garst, B. Binz, F. Jonietz, S. Miihlbauer, 
C. Pfleiderer, and A. Rosch, Phys. Rev. B 86, 054432 
( 2012 ). 

28 L. Kong and J. Zang, Phys. Rev. Lett. Ill, 067203 (2013). 

29 M. Mochizuki, X.Z. Yu, S. Seki, N. Kanazawa, W. 
Koshibae, ,J. Zang, M. Mostovoy, Y. Tokura, and N. Na¬ 
gaosa, Nature Mater. 13, 241 (2014). 

30 A.A. Kovalev, Phys. Rev. B 89, 241101(R) (2014). 

31 S.-Z. Lin, C.D. Batista, C. Reichhardt, and A. Saxena, 
Phys. Rev. Lett. 112, 187203 (2014). 

32 Y.-H. Liu, Y.-Q. Li, and J.H. Han, Phys. Rev. B 87, 100402 
(2013). 

33 J.S. White, K. Prsa, P. Huang, A.A. Omrani, I. Zivkovic, 
M. Bartkowiak, H. Berger, A. Magrez, J.L. Gavilano, G. 
Nagy, J. Zang, and H.M. Rpnnow, Phys. Rev. Lett. 113, 
107203 (2014). 

34 A. Fert, V. Cros, and J. Sampaio, Nature Nanotech. 8, 152 
(2013). 

35 J. Sampaio, V. Cros, S. Rohart, A. Thiaville, and A. Fert, 
Nature Nanotech. 8, 839 (2013). 

36 Y. Zhou and M. Ezawa, Nature Commun. 5, 4652 (2014). 

37 S. Bhattacharya and M.J. Higgins, Phys. Rev. Lett. 70, 
2617 (1993). 

38 A.E. Koshelev and V.M. Vinokur, Phys. Rev. Lett. 73, 
3580 (1994); C.J. Olson, C. Reichhardt, and F. Nori, Phys. 



15 


Rev. Lett. 81, 3757 (1998). 

39 H. Fangohr, S.J. Cox, and P.A.J. de Groot, Phys. Rev. B 
64, 064505 (2001). 

40 A.A. Thiele, Phys. Rev. Lett. 30, 230 (1973). 

41 K. Harada, O. Kamimura, H. Kasai, T. Matsuda, A. Tono- 
mura, and V.V. Moshchalkov, Science 274, 1167 (1996). 

42 C. Reichhardt, C. J. Olson, and F. Nori, Phys. Rev. B 57, 
7937 (1998). 

43 J. Gutierrez, A.V. Silhanek, J. Van de Vondel, W. Gillijns, 
and V.V. Moshchalkov, Phys. Rev. B 80, 140514 (2009); S. 
Avci, Z.L. Xiao, J. Hua, A. Imre, R. Divan, J. Pearson, U. 
Welp, W.K. Kwok, and G.W. Crabtree, Appl. Phys. Lett. 
97, 042511 (2010). 

44 T. Bohlein, J. Mikhael, and C. Bechinger, Nature Mater. 
11, 126 (2012). 

45 C. Reichhardt and F. Nori, Phys. Rev. Lett. 82, 414 
(1999). 

46 P.T. Korda, M.B. Taylor, and D.G. Grier, Phys. Rev. Lett. 
89, 128301 (2002). 

47 M.P. MacDonald, G.C. Spalding, and K. Dholakia, Nature 
(London) 426, 421 (2003). 

48 A. Gopinathan and D.G. Grier, Phys. Rev. Lett. 92, 
130602 (2004). 

49 C. Reichhardt and C. Olson Reichhardt, Phys. Rev. E 69, 
041405 (2004). 


50 A.M. Lacasta, J.M. Sancho, A.H. Romero, and K. Linden- 
berg, Pliys. Rev. Lett. 94, 160601 (2005). 

51 M. Balvin, E. Sohn, T. Iracki, G. Drazer, and J. Frechette, 
Pliys. Rev. Lett. 103, 078301 (2009). S.R. Risbud and G. 
Drazer, Phys. Rev. E 90, 012302 (2014). 

52 C. Reichhardt and C. J. Olson Reichhardt, Phys. Rev. 
Lett. 106, 060603 (2011); T. Bohlein and C. Bechinger, 
Phys. Rev. Lett. 109, 058301 (2012). 

53 H.G. Schuster, Deterministic Chaos (VCH Verlag, Wein- 
heim, 1988). 

54 P. Bak, Phys. Today 39(12), 38 (1986). 

55 S. Shapiro, Pliys. Rev. Lett 11, 80, (1963); A. Barone and 
G. Paterno, Physics and Applications of the Josephson Ef¬ 
fect (John Wiley, New York, 1982). 

56 C. Reichhardt, R.T. Scalettar, G.T. Zimanyi, and N. 
Grpnbech-Jensen, Phys. Rev. B 61, R11914(R) (2000). 

57 E. Scholl, Nonequilibrium Phase Transitions in Semicon¬ 
ductors (Springer-Verlag, Berlin, 1987). 

58 A. Vanossi, N. Manini, and E. Tosatti, Proc. Natl. Acad. 
Sci. (USA) 109, 16429 (2012); D. McDermott, J. Amelang, 
C.J. Olson Reichhardt, and C. Reichhardt, Phys. Rev. E 
88, 062301 (2013). 

59 L. Berger, Phys. Rev. B 2 4559, (1970). 



